Discriminating spin through quantum interference 
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Many of the proposed solutions to the hierarchy and naturalness problems postulate new 'partner' fields to 
the standard model particles. Determining the spins of these new particles will be critical in distinguishing 
among the various possible SM extensions, yet proposed methods rely on the underlying models. We propose a 
new model-independent method for spin measurements which takes advantage of quantum interference among 
helicity states. We demonstrate that this method will be able to discriminate scalar particles from higher spin 
states at the ILC, and discuss application to higher spins and possible uses at the LHC. 

PACS numbers: 



I. INTRODUCTION 

Within the year, the Large Hadron Collider is expected to 
be up and running, granting us at long last access to the scale 
of electroweak symmetry breaking and beyond. One of the 
major puzzles we hope the LHC may provide answers to is 
the hierarchy problem | ll |T| [T| f4 1 : the origin and stability 
of the orders of magnitude gulf between the Higgs vev at 
- 300 GeV and the Planck scale at - lO^^ GeV. Without 
experimental results, theorists over the years have collected 
an impressive aiTay of possible solutions to this problem. Ar- 
guably, the leading contender is supersymmetry Q, but there 
are many others: extra dimensions l6lli7Jli8lll91 lfT0llfTT1lfT2l . 
technicolor |13||2||3|, and little Higgs lfT4l to name a few. 
Many of these models also provide a long-lived, weakly cou- 
pled particle suitable to be a dark matter candidate. 

In many of these possibilities the immediate experimental 
signature from the LHC would be the presence of beyond the 
Standard Model (SM) particles partnered with some or all of 
the known particles. For example, the minimal supersymmet- 
ric standard model (MSSM) doubles the number of particles 
by adding a new field with the same gauge quantum numbers 
and Yukawa couplings as in the SM, but spins differing by 
one half of a unit. Alternatively, the minimal universal extra 
dimensions (UED) model |15| has compactified extra dimen- 
sions which solve the hierarchy problem by 'ending physics' 
at the scale of electroweak symmetry breaking. That is, the 
Planck scale of the true 4 + d dimensional theory is not far 
above a TeV, but appears much larger in 4D once the com- 
pactification occurs. This results in a tower of KK states, each 
containing a heavier version of the SM particles with identical 
quantum numbers, including spin. 

It is well known that, due to the similarities in the particle 
spectrum and quantum numbers, it may be difficult to distin- 
guish the signatures of MSSM particles from the KK=1 modes 
of UED at future collider experiments 1 1 6 1 1 1 7 ] . The existence 
of the KK=2 modes could serve as a discriminator between su- 
persymmetry and extra dimensions, but their high masses may 
make them kinematically inaccessible. Even if produced, they 
typically decay through KK=1 states, and so their presence 
would only be felt through an increase in the KK=1 produc- 



tion cross section ifTSllfTOJI . Determining the spin of the new 
particles will be necessary to confirm the theory underlying 
any new particles. 

There have been several proposals for measuring spin in fu- 
ture collider experiments. The possibilities at a linear collider 
are far more numerous, due to the control over the center of 
mass energy in each event. Threshold scans can distinguish 
scalars from spinors or vector bosons, as the former cross sec- 
tion rises like (3^ while the latter two are proportional to /3 
||20| . However, such a method cannot be used at a hadron col- 
lider, and cannot discriminate between spin 1/2 and spin 1. 
The differential cross section with respect to production angle 
in s-channel pair produced scalars is proportional to sin^ 9, 
while for spinors it is 1 + ^27^2 cos^ 9. Model dependence 
may be present in the form of i-channel diagrams, which in- 
troduce a forward peak which is similar for both spin statistics 
EOl . Such diagrams make the production angle measurement 
of spin more difficult, but may be possible in some cases f2T|. 

The polar angular dependence in decays can also be used 
for spin measurements. However, extracting spin from these 
measurements assumes knowledge of the final state spins and 
also requires chiral couplings, introducing a model depen- 
dence on the spin measurement |22|. While this method was 
originally proposed for the ILC, it was shown that, with suf- 
ficiently long decay chains and exploitation of the asymmetry 
in production of squarks versus antisquarks, supersymmetric 
spinors could be distinguished from phase space decays at the 
LHC 1 1 8 1 1 22 1 1 23 1 1 24 1 1 25 1 1 26 1 . Yet this method reUes heavily 
on the underlying models as the entire decay chain must be 
considered. 

Clearly, determination of spin is a problem still requir- 
ing novel solutions. In this paper we investigate a model- 
independent method to determine the spin of new particles at 
the ILC, first proposed in |[27l . Through interference between 
the different helicity states in a coherent sum, the cross section 
of pair produced particles decaying to two body final states de- 
velops a non-trivial dependence on the azimuthal angle (p of 
the decay. By extracting this dependence, one can determine 
which helicity states entered into the sum, and thus the spin 
of the decaying particle. This method is similar to (and was 
inspired by) the determination of the quark spin measurement 



at SPEAR fIS). At the end of this paper, we will discuss how 
this general method may be extended to the LHC. 

The paper is organized as follows. In section III] we derive 
the angular dependence of the cross section as a function of 
particle spin. We then determine appropriate experimental 
quantities and develop the necessary measurement techniques. 



In section III we apply our method to distinguish scalars in 
pair production at the ILC from production of higher spin 
states. Spin 1/2 and 1 measurements are considered in sec- 
tion|lV]and we conclude in section[Vl Additional calculations 
are supplied in the Appendices. 



II. AZIMUTHAL ANGULAR DEPENDENCE 

To determine the azimuthal dependence of the cross section 
for pair production followed by decay, we start with a particle 
of helicity h moving in the z direction. When this decays into 
a two-body final state, the momenta of the daughter particles 
are confined to a decay plane. If we consider the rotation of 
this plane about the z axis by an angle <j), it is clear that the 
action of the rotation on matrix element of the decay must be 
equivalent to the action of the rotation on the parent particle. 

Rotations of the particle about the z axis introduce a phase 
g-ijj0^ where Jz is the total angular momentum in the z di- 
rection. However, as the momentum is defined to be in the z 
direction, the orbital component is zero, and J^ reduces to h: 



' = m = ^ = 



(1) 



Therefore, the dependence of the decay matrix element 
M decay on (f) must be 

MdecayiK </)) = c'''^ M decay {h, = 0). (2) 

Were we to produce particles in only one helicity state, 
then the total cross section (proportional to the square of 
■Mprod.-Mdecay) would be independent of (p. However, if 
more than one helicity states are produced and then decay, the 
total cross section is proportional to the coherent sum squared: 
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This expression is true only within the validity of the 
narrow-width approximation. However, for 'weakly-coupled' 
physics, such an approximation is justified ||3TI . 

As a result of this interference among the various helic- 
ity states, the cross section develops a cos{n(f)) dependence, 
where n is an integer running from zero to twice the largest 
value of h for the decaying particle. That is to say, the 4> de- 
pendences for a scalar, spinor, and massive vector boson can 
be written as 



cr(s = 0) = 


= ^0 


(4) 
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= Ao + Ai cos(0) 


(5) 


a{s=l) -- 


= Ao + Ai cos(0) + A2 cos(2(/)) 


(6) 



where the Ai are not functions of (/> (though they are non- 
trivial functions of the other kinematics of the problem). The 
exact forms of the cross section must be worked out from the 
standard rules of constructing matrix elements from Feynman 
diagrams, in which case the (j> dependence will become ap- 
parent. However, from this general argument the relationship 
between spin and (j) dependence is made clear 

To exploit this dependence, we consider pair production 
of particles from e+e^ at the future ILC. As motivated by so- 
lutions to the dark matter problem, we expect the production 
of beyond the SM particles to cascade down to some weakly 
coupled particle which will escape the detector Such WIMPs 
are present in the supersymmetric spectrum as the lightest su- 
persymmetric particle (LSP), typically the lightest neutralino; 
or in universal extra dimensions as the lightest Kaluza-Klein 
odd particle (LKP), typically the Bi. Examples of such event 
topologies in the UED and SUSY models are shown in Fig. [T] 
However, our methods do not rely on such specific models. 




FIG. 1 : a) Pair production of KK= 1 muons in universal extra di- 
mensions decaying to opposite sign muons and missing energy in 
the form of two Bi gauge bosons (the LKP). b) Pair production of 
smuons in supersymmetry decaying to opposite sign muons and the 
lightest neutralinos as LSP missing energy 

Measuring the azimuthal dependence of the cross section 
requires that we are able to reconstruct the momentum of the 
parent particle. For simplicity, we specialize to cases where 
the pair-produced particles each decay to a charged lepton and 
missing energy, in which case the events of interest consist of 
£^£^^. While we risk losing some model independence at 
this stage, such signatures are fairly generic in many exten- 
sions to the SM. 

Let the pair produced particles who's spin we wish to mea- 
sure (/ii in Fig. [it or fi in Fig. [TTd) have 4-momenta p^ and 
Pb and mass AI . These particles each decay to visible (effec- 
tively massless) leptons and a weakly coupled particle with 
mass m (the Xi or Bi in Fig. [lb. We label the visible lepton 
momentum pi and p2, and assume that the particles running 
in each leg are identical. The production angle 6 and decay 
angles 6i and (f>i (i — I for the decay of A and i = 2 for the 
decay of B) defined relative to the production plane are shown 
inFig.|2] 

At the ILC, assuming knowledge of the masses M and m, 
it is possible to completely reconstruct the 4-momenta pA and 
Pb (and thus the angles (j>i and (1)2) up to a two-fold ambiguity 
1 32 1 1 33 1 . We note that there are 4 unknown values for both of 
the missing particles in the decay, for a total of 8 unknowns. 
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FIG. 2: The pair produced /i or /ii in tlie lab frame. Tiie beam axis is 
defined as the z axis, with the production angle Q in the x — z plane. 
The z axis is defined to point along the production axis. The decay 
angle 01 is invariant to boosts along z, and so may be defined in 
either the lab frame or the frame of the decaying particle. The angle 
Q\ is defined in the rest frame of fi" / fj,^ . Decay angles 62 and (j)2 
(not shown) are defined equivalently for the fjC^ //i^ . 



There are 4 measured values of the total missing 4-momentum 
p; and for each massive particle in the diagram there is a mass 
relation, for a total of 4 constraints. Therefore, one would ex- 
pect this event to be completely reconstructible. When solv- 
ing the mass relations however, one finds an ambiguity in sign 
when taking a square root, leading to the two-fold ambiguity 
in the reconstructed momentum. For the details of the recon- 
struction, see Appendix [A| 

With less than perfect knowledge of the masses, muon 
momenta and center of mass energy (from beamstrahlung 
||34| and initial state radiation), the true solution will not 
be reconstructed perfectly. At the ILC, masses of lepton 
and gaugino partners are expected to be measured to one 
part per mille Il35lll36l . the tracking resolution as good as 
- 5 X IQ-^{pt/G€W) Wl, and beamstrahlung/ISR should 
be a few percent ll38l . As such, we expect that the errors in- 
troduced in from these effects will be minimal. 

As we have two solutions for the momenta pA and ps, this 
leads to two solutions each for 0i and (/)2. The extracted sig- 
nal in the azimuthal distribution is therefore obtained in the 
combination of the true and false solutions and compare to the 
expected values given in Eqs. (J4|i, (J5]), and (J6]l. In particular, a 
particle of spin n/2 should have Ai — for all i > n. 



III. SCALARS VS. SPINORS 

To demonstrate the utility of this method of spin determi- 
nation, we consider the pair production of scalar light-handed 
smuons in supersymmetry which decay to muons and LSP 



-r.o 



+XiXi)- We compare the az- 
^2 in this scenario to that in 



imuthal distributions of 01 and 
the pair production of ^irs decaying to muons and LKP BiS 
in a UED model (e~e+ — > fJ-i^JJ-iR — * 1^^ IJ-^BiBi). The 
Feynman diagrams for these processes are shown in Fig. [T] 
Analytic formulae for the production and decay cross sections 
for both models are presented in Appendix [B| We stress that 
SUSY and UED are chosen only as benchmark models with 
differing spins and similar final states, the method used to de- 
termine spin can in principle work equally well for any other 





SPS3 


MUED 


X?/Bi 


161 GeV 


302 GeV 


iR/ilR 


181 GeV 


304 GeV 


iL/hL 


289 GeV 


309 GeV 


xt/wt 


306 GeV 


327 GeV 


vtlvii 


276 GeV 


309 GeV 



TABLE I: Relevant particle spectra for the mSUGRA parameter 
point SPS 3 and the minimal universal extra dimension parameters 
R-^ = 300 GeV, A = 20i?~^ and mn = 120 GeV. The MUED 
spectrum was derived using the MUED package [42| for CalcHEP 
1431 . I here refers to the light charged leptons: electrons or muons. 



scenario. 

Representative spectra are required for both supersymmetry 
and universal extra dimensions. In addition, we wish to avoid 
any possible model-specific effects on the azimuthal distribu- 
tions arising from different choices of spectra for supersym- 
metry and extra dimensions. Therefore, as the masses of the 
/i//ii and X1/-B1 are assumed to be known, we perform our 
analysis twice for each model. In the first case we assign the 
masses of the /i and B partners as per a SUSY spectrum, and 
then repeat the processes with the UED case. 

As a representative supersymmetry point, we chose SPS 
3 Il39lll40l in mSUGRA parameter space, which has mo = 
90 GeV, mi/2 = 400 GeV, Aq = 0, tan/3 = 10 and a pos- 
itive [L parameter. Universal extra dimensions is represented 
by the minimal version (MUED) ||4T1| . which requires only 
four parameters: the number and radius R of the extra dimen- 
sions, the scale A to set flavor-universal boundary terms equal 
to zero, and the Higgs mass. We chose one extra dimension 
with i?-i = 300 GeV, A = 20i?-i and a Higgs mass of 
120 GeV. The resulting particle spectra at the TeV scale are 
shown in Table |l] 

Backgrounds consist of the standard model production of 
W~ W^ pair production with leptonic decays to muons and 
neutrinos, ZZ production with decays to muons and neutri- 
nos, and model-background of Xi^xT/^i^^i production 
decaying to muons and i//iyi. While kinematic cuts on the 
invariant mass of the muon pairs can greatly reduce the SM 
background, more efficient cuts can be obtained by requiring 
successful reconstruction of the JIr/ Hir momentum as out- 
lined in Appendix [a] 

The reconstruction algorithm assumes that the masses of 
the produced and escaping particles are known. By assuming 
that the signature of ji^ 11^$ arises from pair production of 
jiR (or Hir) decaying to LSP or LKP, all other events with the 
same signature but different particle masses develop inconsis- 
tencies in their reconstruction. That is, the visible momenta 
are not compatible with the pair production of particles with 
masses other than that for the jir/ iiir, decaying into particles 
with masses other than that of the LSP/LKP. In practice, the 



parameter y defined in Eq. (Alii becomes imaginary. 

With perfect knowledge of masses and muon momentum, 
requiring reconstruction to succeed cuts nearly all of the back- 
ground events. Once detector smearing and mass measure- 
ment errors are included, it is inevitable that some background 
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FIG. 3: Cross sections times brancliing ratios as a function of 
thie beam energy for the UED process e~e"^ —^ MiijMifl ^^ 
fj,~ fj,'^BiBi and tiie SUSY process e~e'*' -^ I^Rp-i ^^ 
/i^/i+XiXi- Figure a) uses the SPS 3 spectrum, while b) uses the 
MUED spectrum (see Table [l]l. 



FIG. 4: Histograms of number of events per azimuthal angle (f> 
for both the true solution to the reconstruction algorithm and the 
combined true and false distribution. The center of mass energy is 
V's — 370 GeV and the luminosity if 500 fb"'^. Figure a) shows the 
UED distribution for e^e^ -^ I^irI^ir ~^ fi~ j-i'^ BiBi while b) is 
the SUSY distribution for e^e^ -^ i^Rf'-R ~^ MM^XiXi- 



will survive the reconstruction cut. Again, with the small er- 
rors in mass measurements available at the ILC we do not 
expect large backgrounds to pollute the data set. 

The total center of mass energy at the ILC is expected to 
reach up to 1 TeV, and an integrated luminosity of 500 fb^ 
is not unrealistic. For the mass spectra chosen, the resulting 
cross sections times branching ratios are shown in Fig. [3] for 
-y/s running from threshold up to 1 TeV. As a result, we expect 
several thousand to several hundreds of thousands of events 
available. 

To simulate the effects of necessary cuts due to the geome- 
try of the detector, we place cuts on the pseudo-rapidity rj. We 
require 77 to be less than 2.5 for both visible muons, as other- 
wise the leptons would vanish unseen down the beam. Also, 
if the missing momentum also points down the beam pipe we 
cannot be sure that the missing energy is truly in the form of 
WIMPs and not merely unobserved SM particles, so we cut 
on T] for missing px as well. 

Using HELAS 1291 the production and decay matrix ele- 
ments were calculated at tree level for each helicity state. Us- 
ing the narrow width approximation, the cross sections as a 
function of 9, 0i, di, (j>2, and 62 were obtained. BASES ll30l . 
an adaptive Monte Carlo, program was used to integrate over 
the other kinematic angles to determine the differential cross 
sections with respect to cf)i and (f>2- As both decaying particles 
have the same spin statistics, the differential distributions are 
the same for both (f)s and so, to increase statistics, the distri- 
butions for 01 and 02 were added. 

Representative distributions for scalar and spinors (includ- 
ing rapidity cuts) are shown in Fig.|4] As can be clearly seen in 
Fig. pk, both the true and false UED distributions have clear 
cos (j) dependence, as expected from spinor decay (Eq. J5])). 
The true distributions for the scalar SUSY decay in Fig.|^is 
flat, as expected from Eq. (J4]). It is therefore apparent even at 
this level of analysis that the (p dependence of the distribution 
contains the spin information necessary for our method. 

Considering the combined true and false distribution in 
Fig. HI a systematic issue for our method becomes readily ap- 
parent. An unexpected cos 20 dependence develops due to 
the false distribution and rapidity cuts, a situation we regard 
as an indication of practical limitations to our method, not 
a fundamental flaw. Whereas the cos 20 dependence may be 
unimportant for the discrimination of scalar versus higher spin 



states, it will become important in distinguishing spinor from 



vectors (section IV 1. Though harder to see by eye, the UED 
distribution also develops a cos 20 dependence in the false so- 
lution. As such, we fit not to Aq + Ai cos but rather to 
Aq + A\ cos + y42 cos 20. 

The overall scaling of the Ai parameters in Eqs. (HI), (Bll, and 
(|6]l depends on the total number of events, which is a function 
of the total cross section. To remove this model dependent 
effect, the parameter of interest in spin determination is not 
Ax, but Ai/Aq, which is independent of the scaling due to 
total cross section. 

Using the least squares method the generated distributions 
were fit to Aq + Ai cos + A2 cos 20. The errors for each pa- 
rameter Ai were obtained after marginalizing over the other 
two parameters. The ratio Ai/Aq for the scalar jiji and spinor 
jj^ijf are shown in Fig. l5] As can be seen, for both the SPS3 
and MUED spectra the values of Ai/Aq for flu are consistent 
with zero for all energies and for both the true and false distri- 
butions. For the spinor /ii/;, the ratio is manifestly non-zero, 
allowing us to distinguish scalars from higher spin states. 

Several aspects of Fig. |5] require closer examination. The 
large error bars for the supersymmetric particles in both spec- 
tra are due to the relatively poor statistics compared to the 
pair production of the spinor KK modes in universal extra di- 
mensions. This is especially apparent near threshold. For the 
spinor particles we also note that, near threshold, the signal 
is on the order of 10%, and decreases towards zero at pro- 
gressively higher energies. This decrease can be readily ex- 
plained as follows: far from threshold, the mass of the pair 
produced particles becomes less relevant, and so their spins 
become more correlated due to chirality conservation. To de- 
termine the distribution of 0i (02), we integrate over all other 
angles in the problem, including 02 (0i). Due to the corre- 
lation of spins in this energy regime, this integration results 
in decoherence of the sum of matrix elements. That is, rather 
than considering \J2h-^i^)\ ' ^^ ^S^ energies the cross sec- 
tion becomes proportional to J^h 1-^(^)1 ' which has no az- 
imuthal angle dependence due to the lack of interference be- 
tween terms. 

Finally, in considering the distribution of true solutions ver- 
sus that of the combined solutions, we note that for the spinor 
case the signal is less once the false solutions are added in. 
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FIG. 5: Top: Ratio Ai/Aq for mSUGRA parameter point SPS3 as 
a function of energy for botii scalar (SUSY) and spinor (UED) pair 
production witii 500 fb~^ luminosity.. Error bars correspond to 95% 
exclusion region. Blue lines correspond to true solution only with 
no rapidity cuts, black dashed lines are true solutions with rapidity 
cuts, red lines to true and false solutions without cuts, and green are 
true and false solutions with cuts. Bottom: Ratio Ai/Ao for MUED 
parameters as in Table ll] for both scalar (SUSY) and spinor (UED) 
production. Color labeling identical to the above. 



FIG. 6: Top: Ratio A2/A0 for mSUGRA parameter point SPS3 as a 
function of energy for spinor (UED) pair production with 500 fb^^ 
luminosity. Error bars correspond to 95% exclusion region. Blue 
lines correspond to true solution only with no rapidity cuts, black 
dashed lines are true solutions with rapidity cuts, red lines to true and 
false solutions without cuts, and green are true and false solutions 
with cuts. Bottom: Ratio A2/A0 for SPS3 parameters for spinor 
(UED) production after correcting for effects of false distribution and 
cuts on a flat distribution. Color labeling identical to the above. 



At low energies the difference between the two is compara- 
tively small, but grows as we move away from threshold. This 
agrees well with the naive intuition that the false distribution 
should be flat in (pi and ^2; however we stress that at higher 
frequencies such intuition fails us and the flat distribution may 
develop non-trivial cos 2(f) dependences. 

To demonstrate this effect we plot in Fig.l6]the ratio A2/A0 
for the decay of spinor /ii/j (using SPS3 parameters). As can 
be seen in the top plot, the true solution without cuts has a co- 
efficient of zero for the cos 2(j> term, as predicted by Eq. (|5]l for 
spinor decay. However once cuts and the false solutions are 
added a non-zero value is generated. Clearly, this can cause 
confusion between a spin- 1/2 particle and a vector or higher 
spin state. 

To attempt to correct for this effect we generate events in 
which the particles decay according to phase space. This flat 
distribution is reconstructed using the method outlined in Ap- 
pendix [a] and rapidity cuts are applied just as in the SUSY 
and UED cases. As a result, the flat distributions also develop 
a cos 20 dependence. The resulting values for A2/AQ using 



only true solutions (with and without cuts) and then both so- 
lutions (with and without cuts) are subtracted from the appro- 
priate spinor ratios to isolate the spin-dependent effect. The 
resulting corrected A2/AQ values are displayed in the lower 
plot in Fig.l6] As can be seen, the flat distribution corrects the 
cos 20 contribution due to cuts but does not remove the false 
distribution's effect, leaving a ^ 0.5% spurious signal at high 
energies. For reasons we do not yet fully appreciate, at low 
energies the false distribution's effects are minimal, allowing 
for the possibility of accurate spin measurements. However, 
it is exactly in this regime that statistics are poor due to the 
proximity of the threshold. 



IV. SPINOR VS. VECTOR 

Due to the large Ai/Aq signal for non-scalars (on the order 
of 10%) and minimal effect of rapidity cuts and false distribu- 
tions on this ratio, the ILC should have little difficulty discern- 
ing that a particle is spin-0. However for higher spins the cuts 



and false solutions introduce potentially dangerous higher fre- 
quency contributions, as has been demonstrated. 

As a result, the question still remains whether this method 
may be practically applied to discriminate spinors from vec- 
tors in general cases. We therefore consider a case of 
pair production of massive vector bosons in UED contrasted 
with spinor production in SUSY. In particular, we consider 
e^ej — > WiWi -^ i^ t'^ Diivui in universal extra dimen- 
sions and e^ej -^ Xi^X]^ ^ l^l'^DpDii in supersymmetry 
where the leptons £ and £' can be either muon or electron type 
(see Fig. It]). Even though the i^i/i' are not the LSP/LKP, they 
decay to neutrinos and the LSP, neither of which is visible in 
the detector. 
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FIG. 8: Cross sections times branching ratios as a function of beam 
energy for the UED process e£ e J -^ W^ W^ 



the SUSY process e^ e]- 



XiXi 
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'/ uii . Figure a) uses 



the SPS 3 spectrum, while b) uses the MUED spectrum (see Tablelll. 




FIG. 7: a) s-channel and b) f-channel pair production of KK= 1 
W bosons in universal extra dimensions decaying to opposite sign 
leptons and missing energy in the form of two i/is. c) s-channel 
and d) f -channel pair production of charginos Xi in supersymmetry 
decaying to opposite sign leptons and sneutrino missing energy 



For these final states, the total cross sections times branch- 
ing ratios as a function of energy are shown in Fig. l8] Once 
again, the supersymmetric cross section is considerably less 
than that in extra dimensions. Furthermore, the small mass 
splittings in the MUED spectrum lead to small cross sec- 
tions compared to the SPS 3 case. Backgrounds include 
SM W^ and ZZ production, and model backgrounds from 
X2X2/^f ^f ' ^"d i^i^/liii pair production decaying to 
charged leptons and missing energy. However, we once again 
find that demanding successful reconstruction effectively cuts 
the background to negligible levels. In addition, we apply the 
T] < 2.5 cuts on the charged leptons and missing momentum. 

We perform fits to Aq + Ai cos <j> + A2 cos 2(j) as in sec- 



tion 



III 



and consider the ratio ^42/^0, using 1 ab~ of in- 
tegrated luminosity (due to the smaller cross sections). The 
results for the SPS 3 spectrum are displayed in Fig. l9j and 
those of the MUED spectrum are shown in Fig. [TO] Note 



that the true solutions for the vector bosons consist of an ap- 
proximately 1% signal in the SPS3 spectrum and ^ 0.5% in 
MUED. In both spectra the true solution for spinors is consis- 
tent with zero. As with the production of fiifi however, the 
presence of the false distribution introduces significant spuri- 
ous values of A2/A0, dwarfing the true signal by a factor of 
-5. 

In the SPS3 spectrum, even with 1 ab^^ the error bars on 
the true solution for the vector bosons barely exclude zero at 
95% confidence. For the MUED case, the situation is much 
worse, as a smaller signal is combined with cross sections sup- 
pressed by nearly an order of magnitude compared to those in 
the SPS3 case. Thus, statistics may be a limiting factor in 
measuring non-zero spins. 

We attempt to correct for the effects of cuts and false so- 
lutions by generating events which decay according to phase 
space. As the production angle may be measured up to the 
two-fold reconstruction ambiguity, we generate the particles 
with the correct 9 distributions and flat 9i, 4>i distributions and 
run the resulting events through the reconstruction and detec- 
tor simulator. The resulting values for A2/AQ are subtracted 
from those in Figs. l9]and 10 in an attempt to isolate the spin 
effects arising from quantum interference from the non-zero 
A2/A0 coming from cuts and the false solutions. The adjusted 
results are shown in Fig. 1 1 for the SPS3 spectrum. Due to the 
small signal and poor statistics in the MUED spectrum, even 
the uncorrected signal in the true solution cannot be distin- 
guished from zero, so we do not adjust for cuts or the false 
solutions. 



Examining Fig. 11 we find that the flat distribution captures 



the effects of cuts on the ratio A2/A0 but does not correctly 
account for the false distributions. We do find that the false 
distributions do not contribute significantly to the ratio near 
threshold, as in the measurements of Ai/Aq. Once again, this 
behavior is not well understood and statistics in this regime 
are limited. It is conceivable that better results would be ob- 
tained by coupling a flat decay in 4^ with the measured dis- 
tribution of 9i to attempt to account for the false distribution. 
This matching has not been performed as yet. 

Thus, while the quantum interference measurement for 
spin-0 stands on solid ground, the situation for higher spins 
is less certain. Even neglecting the issue of false solutions, 
the vector boson ratio A2/A0 is on the order of 1%, and so re- 
quires significant statistics in order to distinguish from spinor 
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vector boson production e£ ej -^ W^W^ -^ £~l''^uieuiii for the 
same spectrum. Color coding as in Fig. [5] Error bars correspond to 
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decays. Furthermore, the false distribution introduces a spuri- 
ous A2/AQ value which has not been fully understood by the 
authors and can dwarf the signal. Finally, the case of of the 
MUED spectrum demonstrates that, while the method of spin 
measurement is model independent, it is vulnerable to model- 
dependent effects such as total cross section, which control 
the statistical error of the fit. However, note that we could 
do much better statistically by adding hadronic final states for 
one of the Xi /W^i while requiring leptonic final states for the 
other. We again would have two-fold ambiguity, but the rest 
of the measurement remains the same as long as we can mea- 
sure the hadronic energies well enough. This may be possible 
by using the energy flow method that matches tracking and 
calorimeter information. 



V. CONCLUSIONS 

We have demonstrated that the quantum interference of 
multiple helicity states can provide a model independent 
method of spin measurements at the ILC. Specifically, with 
reasonable luminosities, scalar particles can be easily distin- 
guished from spin- 1/2 or higher possibilities in pair produc- 



tion followed by decays to visible leptons and missing energy. 
Determining whether a particle is spin- 1/2 or spin-1 suffers 
from two major problems: the first is simply statistics: as the 
signal is on the order of 1%, the requisite luminosity will be 
a stretch for the ILC, at least in the SUSY and UED models 
considered. 

The second issue concerns the false solution to the recon- 
struction of the pair-produced particles' 4-momentum, and 
hence the derived values of the azimuthal angles 4>i and 02 ■ 
With 8 missing momentum components from the two weakly 
interacting particles escaping the detector, 4 measured total 
missing momenta, and 4 mass constraints the system can be 
solved only up to a two-fold ambiguity. While the cos </> dis- 
tribution is flat in the false solution, non-trivial dependences 
on cos 2(j) develop. From explicit calculations, these depen- 
dences appear to be different for flat, spinor, and vector boson 
distributions, and so cannot be subtracted from the combined 
solutions without losing the desired model-independence. 

It therefore behooves us to consider methods for full re- 
construction of the event. If the decay proceeds by emitting 
several visible particles in a cascade of particles with known 
mass down to the LSP/LKP, then we may over-constrain the 
decay, allowing for full reconstruction. In particular, if the 
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pair-produced particles decay to the LSP through an interme- 
diate state, then there would be 6 mass constraints on the sys- 
tem. With only 8 unknown quantities and 4 measured values, 
the false solution is no longer present. Unfortunately, all such 
decays considered by the authors so far have too low a cross 
section to provide useful spin measurements. 

However, such lengthy decay chains raises the possibility 
of applying this method to the LHC. At a hadron collider the 
center of mass energy and frame of reference are unknown 
for a particular parton-parton level event. Thus, only 2 mea- 
sured quantities may be obtained in the event: the components 
of missing transverse momentum f . With a multi-step de- 
cay we obtain 6 mass constraints, combining these with the 
measured |> we can solve the system of 8 missing momen- 
tum components up to the two-fold ambiguity. Additionally, 
the reconstruction algorithm can be used in a modified form 
(33 1 to measure the masses in the decay chain as a neces- 
sary preliminary step to determining the azimuthal angles. As 
the cross section for producing TeV-scale particles with color 
charge at the LHC is very large {e.g. ~ 1 pb for g or q pair 
production |44|) it seems likely that we may obtain enough 
statistics in such a case to at least measure the spin of scalar 



particles if not those of spin 1 or 1/2. This may possibly allow 
discrimination between the gluino and the KK gluon. 

Note that the method we proposed can be used and tested 
already in the Tevatron top quark sample. The interference 
between the two helicity states of the top quark should give 
rise to cos dependence. For the lepton+jet mode, one can 
fully reconstruct the event without a two-fold ambiguity, but 
this suffers from VF+multi-jet background and the not-stellar 
jet energy resolution. The purely leptonic mode has two-fold 
ambiguity but less background and better momentum resolu- 
tion. Run-II should already have enough statistics to attempt 
the study of azimuthal distributions, giving the first direct ex- 
perimental hint on the spin 1/2 nature of the top quark. 



APPENDIX A: RECONSTRUCTION 

The two charged leptons in the events shown in Fig. [T] have 
momenta pi and p2 respectively. We define the perpendicular 
momentum in the event p± — pi x p2- We refer to the pair 
produced unstable particles as A (for /.i^j^ or fij^) and B ip^pj 
or //J). The missing 4-momentum from the decay of A is 
i) , while p is the missing momentum from the decay of B. 
Both the particles escaping the detector have mass to, which 
is assumed to be known. 

Since the pair produced particles A and B (with mass M) 
are back to back, it suffices to solve forpA, as pA = —ps- The 
final state leptons are effectively massless, so Pi = P2 — 0. 
For the massive particles, we have 



Pa = 


-pI- 


= M^ 




(Al) 


f\- 


-fl- 


= to2 




(A2) 


) decays into 


fl (^2^ ^"'^ Pi 


(P2), 




h - 


= PA 


-Pl 






h = 


= PB 


-P2 




(A3) 



At the ILC, the energy of the beams E is known, and for 
pair production the total energy in the event mu st b e split 
equally, so p\—p% — E. Therefore, using Eqs. (Al 1, (A2i 
and (|A3|l we may define the following variables 



Cl 



C2 



PA ■ Pl = -{m^ - M^ + 2Ep°) 


(A4) 


PA-P2^~\{m^~M^ + 2Epl) 


(A5) 


PA-PA = E^- A'P 


(A6) 


Pt -Pj {hi = 1,2) 


(A7) 



62 = 



We can write the momentum pa as 

PA = hPi + t2P2 + yp±- 
Using this definition in Eq. ( |A4| i and ( |A5| l, we find 



(A8) 



Cl 
C2 



tiGii + t2CL\2 
ilO'12 + ^2^22 



t2 = 



a22Cl - ai2C2 

2 
ailO-22 — ^12 

aiiC2 - ai2Ci 

2 
1^110.22 ~ ai2 



(A9) 
(AlO) 



Finally, using Eqs. ( |A6| l, ( |A9| l and ( |A10 i 

62 = (i?aii +2iit2ai2 + ^3022) + y|RL|^ 

y 



±^ 



(i^flii + 2tit2ai2 + tla22) 



\P±\' 



(All) 



The ± sign in this last equation is the two-fold ambiguity in 
the reconstruction. 



APPENDIX B: AMPLITUDES 

The matrix elements for right-handed smuon pair- 
production from polarized e^e+ beams are 



M{e^e 



L^R 






g(-l/2 + ,s^) 
c2^(s-4m|) 



(Bl) 



M{e^e 



R^L 



I 4?77,? 

AflAfl) = (-*e2)W 1 - — ^ sine/x 



Sh-s 



cl^{s-Aml) 



Here, -^s is the center of mass energy and m^ is the mass of 
right-handed smuon. The angle 6 is defined as in Fig.l2] 

The decaying /2* goes to /i^ and a right-handed xi- We 
make the approximation that the neutralino is primarily bino, 
and so the decay matrix element is 



MilL^ ^ M^X?) = -V2g'^i 



M "^1 



(B2) 



Mak- 



Here g' is the hypercharge gauge coupling g' — — 

ing the narrow width approximation, the cross section for the 

four-body final state is simply the incoherent sum over initial 

hehcities 



da = -^}_^\M{ 



^l/r'^r/l 



jl il+)MiJlj^^ Hj^Xi)^ 



L.R 






helicity. Thus, the production matrix elements are 



M(e^e+-.MrflaK«(T)) = ( 



Miele+^^,-J,a)^,tJ,a)) = ( 



Mie],et -^ t^^jiiDf^UD) = ( 



c^(s-m|) 

ie^)(l + cos6') X 
2 

c^(s-to|) 

,,e2)?I^sin0 

^ ^ s{-l/2 + sl,) 
c^(s-m|) 

~ze2)^sin0 



M(e^e+-.Mrfl(TKi.a)) = ( 



Miele+^^l-J,ir)^^tRm = ( 



M(e^e+-.Mrfl(TKi.(T)) - ( 



1 



Sw^ 



c^(s-to|) 



ie )(1 + cos6') X 
2 



s(-l/2 + s2^) 



''){!- cos9) X 

2 
I ^w^ 

c^(s-to|) 

.-,.2mn, . „ 
ie^) — ^sine* 

^sM/2 + 4^ 

c^(s-m|) 

. o, 2to„, . ^ 
le^) — ^sme* 

2 



c^(s-m|) 



(B4) 



Here, t corresponds to right-handed helicity, while [ is left- 
handed. 

In the rest frame of the decaying jii^, there are two possible 
helicities (| and |) decaying to right-handed muons and three 
possible polarization vectors for the Bi (e\, A = ±1, 0). For 
the decay of the /ijjj. the matrix elements are 

Mifi-j,{^) ^ fi],B,i-l)) = 
-^(A^r«(T)-MflBi(0)) = g'^^Jm^-ml^x 



(B3) 



where T is the total width of the fliB.- Note the lack of depen- 
dence on 01 and 02, in accordance with Eq. Q. 

Pair production for right-handed mu- 1 requires four helicity 
combinations for the /iii^^s. Recall that KK states of the chi- 
ral muons are massive particles, as such they can have either 



e+'^i/2sin^ 



e 



"T-Si 






^k 



e-0i/2 3in^ 



Mifi^^ii) ^ fi^B,{+l)) = V2g',/^ 



e-^-^^/^^os^ 



(B5) 



We see here the dependence on the heUcity of the /ii/j as in 
Eq. (pb. Similar equations hold for the decay of /ij'^, with 
01 -> (f>2 and 6*1 -^ 62. 

The total cross section for the event is the coherent sum 
over nm helicities and the incoherent sum over the helicities 
h of the electrons and polarizations A of the KK photons: 



...^ E 



L,R.\\' 



E-^( 



'^l/r'^r/l 



fiiji{h)fi+j^^{h')) 



hh' 
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Once again, F is the total width of fiifi and there is an implied 
momentum conserving 6 function. 
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